Have we already detected astrophysical symptoms of space-time noncommutativity ? 
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We discuss astrophysical implications of K-Minkowski space-time, in which there appears space- 
time noncommutativity. We first derive a velocity formula for particles based on the motion of 
a wave packet. The result is that a massless particle moves at a constant speed as in the usual 
Minkowski space-time, which implies that an arrival time analysis by 7-rays from Markarian (Mk) 
421 does not exclude space-time noncommutativity. Based on this observation, we analyze reaction 
processes in K-Minkowski space-time which are related to the puzzling detections of extremely high- 
energy cosmic rays above the Greisen-Zatsepin-Kuzmin cutoff and of high-energy (~20 TeV) 7-rays 
from Mk 501. In these analyses, we take into account the ambiguity of the momentum conservation 
law which can not be determined uniquely from a mathematical viewpoint. We find that peculiar 
types of momentum conservation law with some length scale of noncommutativity above a critical 
length scale can explain such puzzling detections. We also obtain stringent constraints on the length 
scale of noncommutativity and the freedom of momentum conservation law. 

. ll.30.-j, 95.85.Pw, 96.40.-z, 98.70.Sa. 
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I. INTRODUCTION AND OUTLINE 

Recently, much attention has been paid to extremely high-energy cosmic rays (EHECRs) which have energies above 
that attained in any experimental apparatus on Earth It has been pointed out that these EHECRs provide an 

^ opportunity to investigate space-time properties on very short length scales or very high energy scales. The most 
striking feature is that some of these detections seem to be inconsistent with existing physics, in which such detections 
would be restricted by the Greisen-Zatsepin-Kuzmin (GZK) cutoff ||. That is, if we consider the interaction between 
EHECRs and CMB photons, particles with energy > 7 x 10 19 eV from distant sources cannot reach the Earth. There 
is also another anomalous phenomenon similar to this. That is detections of 7-rays above ~ 20 TeV from distant 
sources (> 100 Mpc) reported in Refs. [f|J|]. These 7-rays are expected to interact with infrared background (IRBG) 
photons and not to reach the Earth in a LI scenario || . In spite of exhaustive research, near sources which can explain 
such detections has not been found. Though there are many attempts explaining these anomalous phenomena, there 
is no absolute solution at present @. 

This could imply an encounter with new physics. Some authors argue that violation of Lorentz invariance (LI) 
might solve EHECRs above GZK cutoff &T2 j l. LI violation might also explain detections of 7-rays above ~ 20 TeV. 
This possibility has been argued in Refs. |13|-|16f|. 

One of the ways to introduce LI violation is to consider space-time noncommutativity with deformed LI, which has 
received attention in recent years since it naturally arises in the contexts of string/M theories [^7|-p2|. It has also 
been argued that space-time uncertainty which comes from a fundamental string scale may be related to space-time 
noncommutativity |f23f . 

Apart from string/M theories, space-time noncommutativity also arises as a result of deformation quantization p4[ . 
Amelino-Camelia et al. [p5|-p7[ considered an interesting toy model called K-Minkowski space-time where noncommu- 
tativity is introduced as [x*,t] = i\x l , where A is a free length scale and the index i runs over 1,2,3 [p9|-|3l"||. They 
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obtained a severe constraint on A through an arrival time analysis of signals from a 7-ray burst [|25|,|26| . If we accept 
this scenario, there is no room for detectable symptoms such as anomalous threshold to explain EHECRs |^| . 

In these papers, the velocity of particles was evaluated using a group velocity formula in the usual Minkowski 
space-time. Here, we derive a more realistic velocity formula based on the motion of a wave packet in K-Minkowski 
space-time. With this formula, we find that the space-time noncommutativity does not affect the velocity of massless 
particles. Motivated by this observation, we analyze reaction processes which are related to both detections of EHECRs 
beyond the GZK cutoff and of ~20 TeV photons. In particular, we pay attention to the momentum conservation 
law which has some ambiguities in this model. We propose to determine the form of the momentum conservation 
law by deciding whether or not space-time noncommutativity is consistent with observations. In fact, we can exclude 
some forms of momentum conservation. Though our approach is purely kinematical, our result will provide a strong 
motivation to consider realistic model of space-time noncommutativity p8| |. Throughout this paper, we use the units 
in which c = Ti = 1 . 



II. k-MINKOWSKI SPACE-TIME 



We briefly review K-Minkowski space-time. The basic commutation relations are 

[x\t] = i\x\ [x\x J ]=Q, 



(1) 



where the indices i,j run 1, 2, 3. We can define differentiation, integration |29| and Fourier transformation in this space- 
time p0[ . In order to define Fourier transformation consistently, the energy E and the momentum p = (pi,p 2 ,P3) of 
a particle form a non-Abelian group G which can be written in a matrix form as, 



(E,p) := 



( e XE pi p 2 P3 
10 
10 

Vo 001 



(2) 



Thus, if we denote the additive operator in K-Minkowski space-time by + to distinguish from the conventional one, 
we can write as 



(E 1 +E 2 ,p 1 +p 2 ) := (E 1 ,p 1 )(E 2 ,p 2 ) 



(Ei+E 2 ,P! 



Pi) 



Following Ref . |2(| , we describe a plane wave as 



, ip X iEt 

4>(E,p) = e ^ e , 

and place the t generator to the right of x generator, i.e., ip(s,p) e lEt e ip ' x . Then, the property 
is found. We can also define the wave in the reverse direction as 

^(Etf)- 1 '■' 



(3) 
(4) 
(5) 
(6) 



which implies that (E^p)^ 1 is an inversion of (E,p). 

Because of these noncommutative structures, modification of Poincare invariance is required to describe physics in 
a covariant way |31j|. The rotation and boost generators can be written as 



Mi 



-%mni J m 



d 



d 

dpn 

A 



^=P^-[-P 



l-e 



2XE 



2X 



d , d 

h XpiPj — . 

opi apj 



(7) 

(8) 



Using (||), a finite boost transformation for the i — 1 direction can be obtained as ||3S 
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Pi = 



tanh(Am) sinh£ 



A[l — tanh(Am) cosh£] ' 

P2 = P3 = 0, 



E = to + — In 
A 



1 — tanh(Am) 



1 — tanh(Am) cosh£ 



where £ is a boost parameter and we choose p — and E = to, i.e., to is a rest mass of the particle, for £ = 0. 
Because of (H), the dispersion relation is altered as 



p 2 e~ XE = K 2 , 



\- 2 (e XE + e- XE -2) 

where if is a constant with the dimension of mass. If we take a rest frame, this can be expressed as 

\- 2 (e Xm + e~ Xm -2)=K 2 . 



(9) 
(10) 

(11) 



(12) 



(13) 



III. THE VELOCITY FORMULA 

Here, we derive a new velocity formula which is one of the main results of this paper. The velocity of the particle 
in the usual Minkowski space-time is 

If we apply this in K-Minkowski space-time, \v\ = e~ XE is obtained for a massless particle, where we used Eq. (|T^). 
This formula, together with the data on 7-rays associated with Markarian (Mk) 421 in Ref. [R3[ l eads to the constraint 
A| < 10 -33 meter |^5|,^(i[^3| . Since this discussion depends crucially on the form of Eq. jl4j), i.e., on what is the 
velocity, we reexamine the group velocity formula by forming a wave packet in K-Minkowski space-time as a more 
realistic situation. For this purpose, we consider infinitesimal changes AE and Ap in E and p, respectively, as a 
result of adding (AE' , Ap') as 

(E+AE\p+Ap') = (E + AE,p + Ap). (15) 
In this case, we can express (AE' , Ap') as 

(A£',Ap') = (A£,^§), (16) 
where we keep only terms in first order in AE and Ap. By using Eqs. (|) and ©, we make a wave packet as follows 

I = "0(B-AB,p-Ap) + 1p(E+AE,p+Ap) 

— V>(-E,p)(-A.E',-Ap') + 1p(E,p)(AE' } Ap') 

= y>(£,p)V'(-A.E' ,-Ap') + 1p(E,p)tp(AE',Ap') 



, r -iAp' X -iAE't . iAp' X iAE'ti 

V(E,p)[e e +e ^ e J 

Ap / e XE AEt\ 



^2e i P x e mt cos 

By considering \I\ 2 , the group velocity Vi can be written as 

,d 
dp 



Ap J 



(17) 



v r .= e XE — . (18) 



We also consider a similar relation 

(AE'+E,Ap'+p) = (E + AE,p + Ap), (19) 
which is different from Eq. (flq) due to noncommutativity. In this case, the corresponding group velocity v r is 
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dE 

V r := (20) 

Using (|l2|) and (|l3|), we obtain the important conclusion that massless particles move in a constant speed \vi\ — \v r \ = 1 
as in the usual Minkowski space-time for arbitrary A |3{| . Therefore, the argument in Ref. |2j||2(|^3| does not apply. 
In this case, there appears the possibility that the large value of A (> 10 _33 m) may solve the puzzling problems of 
EHECRs above GZK cutoff and of ~20 TeV photons simultaneously. We investigate this possibility next. However, 
we emphasis on the importance of the result not because K-Minkowski space-time can avoid the constraint but because 
our result provides an opportunity to reconsider LI deformation models in general. 



IV. THRESHOLD ANOMALY 



We first consider the two-body head-on collision of particles and subsequent creation of two particles 1 + 2^3 + 4. 
We define the energy Ei and momentum p i of the i-th particle as those in the laboratory frame. We denote the rest 
mass of the i-th particle as m;. We also assume that m 2 = 0, 7713 7^ 0, TO4 7^ and p i = (pi, 0,0). In the usual 
Minkowski space-time, we use the dispersion relation 

Ef-p^ml (21) 

and the energy momentum conservation law, 

Pl+P2=P3+P4, (22) 

E 1 +E 2 =E 3 + E i , (23) 

to obtain the threshold value of E\, which we denote by E t h,o- We assume that the resultant particles are at rest in 
the center-of-mass frame in the threshold reaction. In the laboratory frame, this means that the resultant particles 
move in the same speed, that is 

PL=P±. (24) 
ma TO4 

We also assume that p 2 has an opposite sign against that of p\. If we neglect higher order terms in E 2 , then 

We also define the threshold value of p\ as pth,o which can be approximated as Pthfi ~ E t h,o- 
Next, we consider the same reaction in K-Minkowski space-time. Eq.(|2l|) is replaced by 

\- 2 (e XE > +e- XE * -2)-( Pl ) 2 e- XE > 

= \- 2 (e Xm ' +e- Xm * - 2). (26) 

If we interpret the algebra in K-Minkowski space-time faithfully, the energy momentum conservation law is expressed 
as 



(E uPl )(E 2 ,p 2 ) = (E 3 ,p 3 )(E A ,p 4 ). (27) 

Even if it holds, one should note that we need a rule to distinguish two particles. If we consider the collision of two 
particles with A, B G G, respectively, does it correspond to AB, BA or anything else ? At present, we have no 
way to determine it. Amelino-Camelia et al. ]26| ] proposed to find the rule by experiments. Here, we introduce a 
phcnomenological parameter a, which controls the form of conservation law as follows: 

a(E 1 ,p 1 )(E 2l p 2 ) + (1 - a)(E 2 ,p 2 ){E uPl ) 
= a{E 3l p 3 )(E ilPA ) + (1 - a)(£; 4 ,p 4 )(^3,P3) • (28) 

As regards plausible values for a, care must be taken. If we consider two particles of the same species, a = 1/2 would 
be physically reasonable value, since if they have same energy and move opposite direction each other, they have zero 
total momentum only for this choice. In fact, the parameter a may be a function of physical quantities of two particles 
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such as mass, charge and/or spin for two different species. Here, we use the same value of a on the left and the right 
hand sides of for convenience. Moreover, we restrict our attention to < a < 1 for clarity. 

We also need to impose the condition that the resultant particles are at rest in the center-of-mass frame. To obtain 
a relation between momenta ps and P4, we use the boost transformation (|^). For the same value of £, we obtain 

?1 = 2± . (29) 

tanh(Am3) tanh(Am4) 

We can solve E\ as a function of a, A, mi, 7713, 771,4 and E 2 by using (|26|), ( p8| ) and (|29|). We apply this result to two 
astrophysical cases. 



A. Threshold anomaly for TeV 7-rays 

Here, we consider the process 7 + 7 — > e + + e - , which may occur when a 7-ray travels in the IRBG. In this case, 
rni = and 7713 = ?ri4 = m e , where m e is the electron mass. If we assume the existence of IRBG photons (0.2 < E 2 5 
eV) then the threshold is E t h,o ~ 1 TeV in Minkowski space-time. Then, the reported detection of ~20 TeV photons 
from Mk501 (^150 Mpc from the Earth) would be difficult to explain 0]. 

We summarize the equation for the threshold in K-Minkowski space-time which is derived from (pfl) , (E8T) and (|2£ 



AB = yx(yx + l) 2 sinh 2 , (30) 

where 

A := (1 - a)y 4 - (1 - 2a)y 2 - a, (31) 

B := ax A + (1 - 2a)ir 2 + a - 1, (32) 

and x := e XEl ^ 2 and y := e XE2 ^ 2 . Since we are considering the collision of two particles of the same species, a — 1/2 
would be physically reasonable. Note that, though we have E t h,o ~ Pth.o for high energy particles in the usual 
Minkowski space-time, this is not the case in K-Minkowski space-time. 

We should recall that, to estimate the energy of primary particles, we calculate the sum of energy of all secondary 
particles. Since energy is conserved in the usual sense even in K-Minkowski space-time, we do not need to take into 
account the effect of space-time noncommutativity to estimate the energy of primary particles. Thus, the observation 
of ^20 TeV photons in usual Minkowski space-time has the same meaning also in K-Minkowski space-time. On the 
other hand, the usual sum of momenta of all secondary particles does not coincide with the momentum of the primary 
particle in this space-time. Therefore, if pth could be evaluated independently of the observation of E t h, it might 
become important to extract information about space-time noncommutativity through the detection of violation of 
the usual momentum conservation. We exhibit properties of both the energy and the momentum from this reason. 

We first show the dependence of E t h and pth on A > in Figs, [j] (a) and (b), respectively. For simplicity, £2 is 
chosen as Ei — 1 eV for IRBG photons. For a = 0, E t h and pth increase with A, compared with the same quantities 
in Minkowski space-time. In particular, E t h and pth diverge for A := A c ~ 4 TeV -1 . That is, the universe is entirely 
tranparent for A > A c . For a — 1/2 and 1, pth increases with A, though E t h decreases. 

For all a, a first-order correction in A arises for p t h- If we expand p t h as p t h — J2T=o tt^i tne first-order coefficient 
Pth, 1 is written as 



/<//>.! = : />//,. Pth,o(l - a) + E 2 (a - i 
On the other hand, the first-order correction in A for E t h, which we denote E t h.i, is written as 



(33) 



E th ,i = E tht0 (E th!0 ~ E 2 ) Q - «) • (34) 

Thus, it disappears for a = 1/2. 

The reason why E^ and p t h disappear for a = above A c ~ 4 TeV -1 can be understood as follows. For AEth ^ 1, 
and Am e , XE 2 < 1, we can approximate eq. PC]) as 
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(1 + XE 2 )x w A£ 2 for a = 0, (35) 

aE 2 x PS ^ for a ^ 0. (36) 



In this range of approximation, since XE t h,o — Xm? e /E 2 is expected to be larger than 1, eq. (p5|) has no real solution 
E t h while a real solution Bj/, exists for a ^ 0. This means that A c for a = is characterized by 1/Eth.o ~ 1 TeV -1 . 

For XEth, Am e , A£?2 S> 1, we can summarize the results as follows. For a = 0, cq. (pQ ) is approximated as 
y 4 a; 2 ~ (ya;) 3 e A ™ lc /4, which yields i?i = — 2m e + i?2 < 0. This contradicts the first assumption. So a solution does 
not exist. In a similar way, we can show that E\ approaches 2m e + E 2 and 2m e — E 2 for a = 1 and for a ^ 0,1, 
respectively. 

To investigate properties for A < 0, we replace A with — A. In cq. (|30|), this corresponds to the replacement x — * 1/x 
and y — * 1/y. We find that eq. ( |30| ) becomes invariant if a is also replaced by (1 — a). 

Thus, the case a <C 1, A > 4 TeV" 1 , and the case 1 - a C 1, -A ^ 4 TeV -1 remain as candidate solutions for 
-20 TeV photons. On the other hand, we can exclude a = 0(1) and A > 10 TeV" 1 , or (1 - a) = 0(1) and -A > 10 
TeV -1 from the present experimental data. 



B. Threshold anomaly for GZK cutoff 

Here, we consider the interaction of ultra high energy protons with CMB photons (~ 10~ 3 eV) which results in a 
pair production p + 7 — > p + ttq. In this case, mi = 1TI3 = m p and 7714 = m„, where m p and m v are the proton mass 
and the pion mass, respectively. Because of Ethfi ~ 7x 10 19 eV, it is difficult for EHECRs above E t h,o to reach the 
Earth from cosmologically distant sources. 



We solve (|26j), (28) and (29) numerically. We show the dependence of E t h and p t h for A > 0, in Figs, g (a) and (b), 
respectively. E 2 is chosen as E 2 = 10 -3 eV for CMB photons. Compared with Fig. |l|, we find that the qualitative 
features for small A are quite similar, i.e., E t h/E t h.a > 1 for a = and E t h/E t h.o < 1 for a = 1, while pth/pth.o > 1 
for all cases. 

However, we find qualitative differences from Fig. 0for A > 10~ 8 TeV- 1 . The threshold disappears for A > A c — 
2 x 10~ 8 TeV -1 in the a = case, which can be explained as in the 7-ray case since A c coincides approximately with 
l/E t h,o- For the case a = 1/2 and 1, E t h/E t h,o increases with A(> 3 x 10~ 8 TeV -1 ) and disappears for A > 5 x 10~ 8 
TeV" 1 , unlike the 7-ray case. 

In this case, there is no simple symmetry about A — > — A, as found in the previous case. Thus, we also show the 
dependence of E t h and p t h for A < 0, in Figs. || (a) and (b), respectively. We have a crucial difference from the 7-ray 
case even for A > — 10~ 9 TeV -1 . For a = 1, there is a value E t hi over which the reaction does not occur. We denote 
Eth2 by a dot-dashed line. E t hi diverges as A — > — and merges with E t h at A = A c — — 7 x 10~ 9 TeV -1 . 

We find that the behavior for small |A| is that E t h/E t h,o > 1 for a = 1 and E t h/E t h,o < 1 for a — as in the 7-ray 
case for A < 0. For a = 1/2 and 1, E th decreases with A(< -5 x 10~ 8 TeV -1 ). 

Thus, the a < 1 case for A > 2 x 10~ 8 TeV" 1 , the a = 0(1) case for A > 5 x 10~ 8 TeV" 1 and the (l-a)«l case 
for A < —7 x 10~ 9 TeV^ 1 remains as candidate explanations for detections of super GZK events. For (1 — a) = O(l), 
we can exclude A < — 10~ 7 TeV -1 . 



V. CONCLUSION AND DISCUSSION 



We have first considered a velocity formula to describe the particle motion based on the motion of a wave packet in k- 
Minkowski space-time. In this formula, space-time noncommutativity does not affect the motion of a massless particle. 
Thus, an arrival time analysis of 7-ray bursts in Refs. (2^j2(^,|33| does not exclude space-time noncommutativity in 
this model. Since this feature had not been discussed so far, it should be stressed and is one of our main conclusions 
here. 

Based on this consideration, we have obtained threshold values for reactions 7 + 7 — > e + + e~ and p + 7 — > p + ttq in 
K-Minkowski space-time and analyzed their relavance to the puzzling observations of —20 TeV photons and EHECRs 
above the GZK cutoff, introducing a parameter a to take into account the ambiguity of the momentum conservation 
law. 

In the TeV 7-ray case, though a = 1/2 is favorable in the physical context, only a <C 1 for A > 4 TeV -1 , or 
(1 - a) < 1 for A < —4 TeV -1 appear able to explain the detections of 7-rays above ~20 TeV. The possibilities 
a = 0(1) for A > 10 TeV~\ or (1 - a) = O(l) for A < -10 TeV" 1 , are excluded. 
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In the EHECR case, we cannot assign definite values to a, because it may depend on, e.g., masses and/or charges 
of two particles. The possibilities a < 1 and A > 2 x 1CT 8 TeV -1 , or a = 0(1) and A > 5 x TeV -1 remain 

viable. We can exclude cases in which (1 — a) = O(l) and A < — 10~ 7 TeV -1 . 

Thus, a <C 1 for A > 4 TeV -1 or (1 — a) -C 1 for A < —4 TeV -1 appear able to explain both phenomena. Our results 
are important because they suggest that extremely high-energy particles might be expected in realistic models with 
space-time noncommutativity. If this is the case, then we might have already detected symptoms of the space-time 
noncommutativity. 
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(b) 

FIG. 1. Threshold anomaly for TeV-7 rays for A > 0. (a) \-E t h, (b) \-p t h are plotted for E% = 1 eV. For a = 0.5 and 1, 
Eth decreases with A increases for A > 1 TeV -1 , while p t h monotonically increases. The a <C case is only desirable to explain 
~20 TeV photons. It is noted that E t h is invariant under the transformation A — > —A and a — > (1 — a). 
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FIG. 2. Threshold anomaly for GZK cutoff for A > 0. (a) X-E th , (b) \-p th are plotted for E 2 = 1CT 3 eV. 
features for small A are similar to those of Fig. 1, they show drastic difference from Fig. 1 for A > 3 x 10~ 



Though qualitative 
1 TeV" 1 . 
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FIG. 3. Threshold anomaly for GZK cutoff for A < 0. (a) X-E th , (b) X-p th are plotted for E 2 = 1CT 3 eV and A < 0. Unlike 
the case A > 0, the threshold vanishes only the a — 1 case for A < —7 x 10 -9 TeV -1 . 
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